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Abstract 

We study the scaled-energy (xb) distribution of bottom-flavored hadrons (B) 
inclusively produced in top-quark decays at next-to-leading order (NLO) in the 
general-mass variable-flavor-number scheme endowed with realistic, nonperturba- 
tive fragmentation functions that are obtained through a global fit to e + e _ data 
from CERN LEP1 and SLAC SLC exploiting their universality and scaling viola- 
tions. Specifically, we study the effects of gluon fragmentation and finite bottom- 
quark and i?-hadron masses. We find the NLO corrections to be significant. Gluon 
fragmentation leads to an appreciable reduction in the partial decay width at low 
values oi xb- Hadron masses are responsible for the Iow-xb threshold, while the 
bottom-quark mass is of minor importance. Neglecting the latter, we also study the 
doubly differential distribution d 2 T / {dx b d cos 9) of the partial width of the decay 
t —> bW + —> Bi + vi + X, where 9 is the decay angle of the charged lepton in the 
TV-boson rest frame. 

PACS numbers: 12.38.Bx, 13.85.Ni, 14.40.Nd, 14.65. Ha 
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1 Introduction 



Among other things, the CERN Large Hadron Collider (LHC) is a superlative top factory, 
producing about 90 million top-quark pairs per year of running at design energy v^S* = 
14 TeV and design luminosity C = 10 34 cm _2 s _1 in each of the four experiments pp. 
This will allow us to determine the properties of the top quark, such as its mass mt, total 
decay width T t , branching fractions, and elements Vt q of the Cabibbo-Kobayashi-Maskawa 
(CKM) [2] quark mixing matrix, with unprecedented precision. Due to its large mass, 
the top quark decays so rapidly that it has no time to hadronize and passes on its full 
spin information to its decay products. If it were not for the confinement of color, the 
top quark could, therefore, be considered as a free particle. Due to \V t b\ ~ 1, top quarks 
almost exclusively decay to bottom quarks, via t — > bW + . 

On the other hand, bottom quarks hadronize, via b — > B + X, before they decay, so 
that the decay process t — > BW + + X is of prime importance, and it is an urgent task 
to predict its partial decay width as realistically and reliably as possible. Of particular 
interest are the distribution in the scaled S-hadron energy xb in the top-quark rest frame, 
and, in the case of leptonic J¥-boson decays W + — > £ + ise, the one in the charged-lepton 
decay angle 9 in the W-boson rest frame. In fact, the xb distribution provides direct 
access to the 5-hadron fragmentation functions (FFs), and the cos# distribution allows 
us to analyze the VF-boson polarization and so to further constrain the 5-hadron FFs 
by exploiting xb distributions for all the VT-boson polarization states. B mesons are, for 
instance, cleanly identified by their decays to J ftp mesons, which are easy to tag through 
their spectacular decays to e + e~ and pairs. 

The theoretical aspects of top-quark physics at the LHC are nicely summarized in a 
recent review paper |3j. In the approximation of treating the bottom quark as a stable 
final-state particle that does not hadronize, the QCD corrections to V(t — > bW + ) are 
known at NLO with subsequent W + —> t^vi decay [415] and at next-to- next-to-leading 
order for stable W boson [6]. The terms of order /3o«™ +1 , where /3q is the first coefficient 
of the QCD beta function and a s is the strong-coupling constant, were resummed to all 
orders in Ref. |7j. The NLO electroweak corrections were found in Ref. [8]. 

The hadronization of the bottom quark was considered in the NLO QCD analyses of 
top-quark decay in Refs. [9|10|lllj and was, in fact, identified to be the largest source of 
uncertainty in the determination of the top-quark mass. In Refs. [9|10|lllj . the W boson 
was taken to be stable, the bottom-quark mass was neglected at the parton level, and 
the hadronization process b — > B + X was implemented as a convolution of a pertur- 
bative FF [T2], describing in a way the conversion of the massless bottom quark to its 
massive counterpart, and a nonperturbative FF modeling the hadronization of the latter. 
Hereby, the perturbative FF depends on a factorization scale and is subject to Dokshitzer- 
Gribov-Lipatov-Altarelli-Parisi (DGLAP) [13] evolution, while the nonperturbative FF is 
scale independent. Because of the treatment of a heavy quark as a massless parton, this 
framework corresponds to the zero-mass variable-flavor-number (ZM-VFN) scheme. In 
Refs. PHU], also soft-gluon resummation was studied. In Ref. [TT] . also the distribution 
in the invariant mass mse of the B hadron and the charged lepton £ from ly-boson decay 
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was considered. 

In this paper, we revisit 5-hadron production from top-quark decay working at NLO in 
the general-mass variable-flavor- number (GM-VFN) scheme, which was elaborated for in- 
clusive heavy-flavored-hadron production in e + e~ annihilation [H], two-photon collisions 
[15] . photoproduction [16], and hadroproduction [TT1[T8] . and provides an ideal theoretical 
framework also here. Being manifestly based on Collin's hard-scattering factorization the- 
orem appropriate for massive quarks [19], this factorization scheme allows one to resum 
the large logarithms in nib, to retain the finite-m?, effects, and to preserve the universality 
of the FFs, whose scaling violations remain to be subject to DGLAP evolution. In this 
way, it combines the virtues of the ZM-VFN and fixed-flavor-number (FFN) schemes and, 
at the same time, avoids their flaws. It is, in fact, a tailor-made tool for global analyses of 
experimental data on the inclusive production of heavy-flavored hadrons, allowing one to 
transfer nonperturbative information on the hadronization of quarks and gluons from one 
type of experiment to another and, within one type of experiment, from one energy scale 
fip to another, without the restriction fip 3> m& inherent to the ZM-VFN scheme. In the 
GM-VFN scheme, the perturbative FFs enter the formalism via subtraction terms for the 
hard-scattering cross sections and decay rates, so that the actual FFs are truly nonpertur- 
bative and may be assumed to have some smooth forms that can be determined through 
global data fits. In contrast to the FFN scheme, the GM-VFN scheme also accommodates 
FFs for gluons and light quarks, as in the ZM-VFN scheme. 

Specifically, our analysis is supposed to enhance those of Refs. [9|10JI11] by retaining all 
nonlogarithmic raj, terms of the result in the FFN scheme and by including light-parton 
fragmentation. Furthermore, we include finite-m^ effects, which modify the relations 
between partonic and hadronic variables and reduce the available phase space. Although 
these additional effects are not expected to be truly sizable numerically, except for certain 
corners the phase space, their study is nevertheless mandatory in order to fully exploit 
the enormous statistics of the LHC data to be taken in the long run for a high-precision 
determination of the top-quark properties. We also extend Refs. [9|10|lll] by including 
subsequent leptonic decays of the W boson and studying the distribution in the angle 9 
of the charged lepton in the P^-boson rest frame, while Ref. [TT] is focused on the m B £ 
distribution. 

This paper is organized as follows. In Sec. [2j we explain how to incorporate finite-m^ 
corrections in the evaluation of dT(t — > B + X)/dxs- In Sec. [3j we give the parton-level 
expressions for dT(t — > B + X)/dxs at NLO in the ZM-VFN and FFN schemes and 
combine them to obtain those in the GM-VFN scheme. In Sec. HI we list the parton-level 
formulas needed to evaluate d 2 T(t — > Blv^ + X)/(dxBdcos6) at NLO in the ZM-VFN 
scheme. In Sec. |5l we present our numerical analysis. In Sec. |6l we summarize our 
conclusions. The Appendix accommodates some formulas that are too lengthy to be 
displayed in Sec. HJ 
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2 Hadron mass effects 



We consider the decay process 

t(pt) -> b(p b ) + W + (p w )(+g(p g )) -> B(p B ) + X, (1) 

where X collectively denotes the unobserved final-state particles and the four-momentum 
assignments are indicated in parentheses. The gluon in Eq. (OQ) contributes to the real 
radiation at NLO. Both the b quark and the gluon may hadronize to the B hadron. 

In the top-quark rest frame, the b quark, gluon, and B hadron have energies E{ = 
p t ■ Pi/m t (i = b,g,B), which range from E™ m = m&, E™ n = 0, and E™ n = mg to 
^ma* = ( m 2 + m 2 _ m ^/{2m t )^ E%** = [m 2 - (m b + m w ) 2 ]/{2m t ) } and E™ ax = {m 2 + 
m 2 B —m 2 v )/(2m t ), respectively. In the case of gluon fragmentation, g — >■ B, the maximum 
i?-hadron energy is E B ax = [m 2 + m 2 B — (m& + mvK) 2 ]/(2m t ). It is convenient to introduce 
the scaled energies Xj = Ei/E™* (i = b,g,B). 

We wish to calculate the partial decay width of process (pP) differential in xb, dT/dxB, 
at NLO in the GM-VFN scheme taking into account finite-m^ corrections. In the ZM- 
VFN scheme, the four-momenta of the produced hadron and the mother parton are related 
as pb = zp a (a = b,g), where the scaling variable z takes the values < z < 1. This 
simple relation is not compatible with finite quark and/or hadron masses and needs to 
be generalized when passing from the ZM-VFN scheme to the GM-VFN scheme. There 
is some freedom in defining the scaling variable in the presence of quark and/or hadron 
masses. In the case under consideration, a convenient choice is Eb = zE a [TJ], i.e. to 
retain just one of the four equations ps = zp a . By the factorization theorem of the 
QCD-improved parton model, we then have 



dT = ^ / dz dt a (jJL R ,fi F ) 
. Jo 

a=b,g 



D a {z,fJL F ), (2) 

E a =E B /z 



where dt a (fiR, fi F ) is the differential decay width of the parton-level process t — > a + X, 
with X comprising the W boson and any other parton, and D a (z,fi F ) is the FF of the 
transition a — >■ B. Here, and fi F are the renormalization and factorization scales, 
respectively. Substituting dT = dxB{dx a / dxs)dT / dx a and eliminating z = xs/x a as the 
integration variable, we obtain our master formula 

dT y-^ [ x ™* dx a dt a (x B \ /Q s 

~d~x~ = 2^ \ ~dx~ im,fIF)Da [~^ F r (3) 

UJjB , J x mln UX a \ X a I 

a=b,g a 

Using xb < x a along with the above bounds on Ei (i = b,g,B), we have x™ in = 
max(pb, xb), x™ ax = 1, x™ n = x B , and x™ ax = p, where = mj/i?™ ax (i = b,g,B) and 
p = e^/E™*. The kinematically allowed x B ranges are p B < x B < min(l, Eg Lax /£'™ ax ) 
for a = b and ps < xb < p for a = g. In reality, we have nib < m B , so that x™ m = xb 
and x^ ax = 1 for a = b. 
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In order to assess the theoretical uncertainty due to the freedom in the choice of scaling 
variable in the GM-VFN scheme, we also consider here the definition in terms of the plus 
component V + = (V° + V 3 )/\/2 of a four- vector V in light-cone coordinates. Taking 
the three-axis to point along the common flight direction of parton a and the B hadron, 
we define p B = zp^ [20]. This definition is invariant under boosts along the three-axis. 
Starting from the factorization formula © with E a = Eb/z replaced by p+ = p%/z, we 
obtain 

i f a r dr a/ . ^ . 

dx B y/x 2 B -f? B ^J x ™ dx a 



where 



Xa + y% a Pa 

and it is understood that m g = 0. Using again the above bounds on E{ (i = b, g, B), but 
imposing z < 1 instead of xb < x a , we now have 




P 2 



xb + \ x% -p%+ %=== , (6) 



a^B + a/ 



-v . — 



while x™ ax and x™ ax go unchanged. The kinematically allowed xb ranges now are ps < 
a^B < E^ ax /E^ ax for a = 6 and p B < x B < E^/E™* ioi a = g. 

Clearly, both Eqs. (J3J and (jlj) may also be used for m h = 0, to improve the ZM-VFN 
scheme by accommodating hadron-mass corrections. If also itlb = is put, then Eqs. ([3]) 
and (jlj) coincide reproducing the familiar factorization formula of the massless parton 
model. 

Taking a closer look at Eq. ([3]), we observe that the GM-VFN prediction implemented 
with the energy scaling variable z = Eb/ E a is not affected by finite-m^ corrections inside 
the xb region accessible for > mj. This is quite different for the implementation 
with the light-cone-momentum scaling variable z = p^/Pa y i a Eq. (jlj). Of course, these 
observations carry over to the ZM-VFN case of m& = 0. 



3 Analytic results for dT(t — ► B + X)/dx B 

We now discuss the evaluation of the quantities dT a (ix R) p F )/dx a at NLO in the GM-VFN 
scheme. Their counterparts in the FFN scheme contain the full m& dependence. However, 
in the limit m&/m t — > 0, they develop large logarithmic would-be collinear singularities 
of the type (a s /n) \n(mf/ml), which spoil the convergence of the perturbative expansion. 
There is no conceptual necessity for FFs in this scheme. They may still be introduced by 
hand, but there is no factorization theorem to guarantee their universality. 

In the ZM-VFN scheme, where = is put right from the start, such collinear singu- 
larities are regularized by dimensional regularization in D = 4 — 2e space-time dimensions 
to become single poles in e, which are subtracted at factorization scale p, F and absorbed 
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into the bare FFs according to the modified minimal-subtraction (MS) scheme. This 
renormalizes the FFs, endowing them with fip dependence, and generates in dT a /dx a 
finite terms of the form (a s /7r) ln(m 2 / fi 2 F ) , which are rendered perturbatively small by 
choosing fip = 0(m t ). In this scheme, m;, only sets the initial scale /z™ 1 = 0(rrib) of 
the DGLAP evolution, where ansaetze for the z dependences of the FFs D a (z,[i l p l ) are 
injected. The DGLAP evolution from fi 1 ™ to fip then effectively resums the problematic 
logarithms (a s /ir) \n(m 2 /m 2 ) of the FFN scheme. However, all information on the mj 
dependence of dT a /dx a is wasted. 

The GM-VFN scheme is devised to resum the large logarithms in rrib and to retain the 
entire nonlogarithmic rrib dependence at the same time. This is achieved by introducing 
appropriate subtraction terms in the NLO FFN expressions for dT a /dx a , so that the NLO 
ZM-VFN results are exactly recovered in the limit m^jm t — > 0. These subtraction terms 
are universal and so are the FFs in this scheme, as is guaranteed by Collin's factorization 
theorem [TJJ]. If the same experimental data are fitted in the ZM-VFN and GM-VFN 
schemes, the resulting FFs will be somewhat different. 

For the sake of a comparative analysis of the GM-VFN and ZM-VFN schemes, we 
need to know dT a /dx a at NLO in both schemes. The NLO expressions for dTb/dxb in the 
ZM-VFN and FFN schemes may be found in Eqs. (5) and (6) and Eq. (A.10) of Ref. [9], 
respectively. We verified them by an independent calculation and also derived those for 
dYg/dxg. For the reader's convenience, we list here all our results. In the ZM-VFN 
scheme, we have 



1 dT b 

r dx b 



1 dTr 



where 



0(1 - x b ) H C F 



1 + x 2 



x 



X 



2tt ^ \1 — Xb 

w(l —w)(l — Xb) 2 



ln x 2 (l - w) 2 m 2 i 



/4 



+ 



4x h 



(l + 2w)[l-x b (l-w)\ 



0+2(1 + 2* 



- 5(1 - x b ) 



4Li 2 (w) + 21nwln(l - w) + 
1 + (1 



2w 



In w + 



l-t 



2tt 



1 — w 

(1 — x g ) 2 (l — w) 2 m 



/ ln(l - x b ) 

5 + Aw 

mil — w ) H 

l + 2w v ' 2 
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i 



2(1 + 2w)[l -x g (l - w)) : 



[1 - Xg(l - W)]fl 

-(1 - 6w)(l - wfx] + 2(1 - w)(l - 2w) 

l + 2w 



x (3 + w)x 2 - (13 - 25w 2 + 6w d )x g + 4(1 + 2w)(3 - 2w) - 4 



,(7) 



G F ml\V tb \ 2 
StiV2 



(1 - w) 2 (l + 2w) 



is the total decay width at LO, Gp is Fermi's constant, w = mly/m 2 , Cp = (N 2 — 
1)/(2N C ) = 4/3 for N c = 3 quark colors, and Li 2 (x) = — J Q (dt/t) ln(l — t) is the Spence 

function. Integrating dTb/dxb of Eq. <^ over x b , we recover the familiar result given by 
Eq. (2.8) in connection with Eqs. (3.1) and (3.2) of Ref. [I]. 
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In the FFN scheme, we have 



1 dT b 
T dx b 



8(1 -x b ) + 



C F a s (fi R ) 
ttQ 



2s 



Li 5 



2Q 



l-s + Q 



— Lie 



2Q 



s-b + Q 
s-b + Q 



111(5 + g) I 111 + ln 2^W) ) + 2 ln 6 ln 2,(1 - 

3^ + S -^ln^ + (l- & )lnI- S + Q 



Gn 



+ Q 
- 2$(x b ) 



Vb 

Q (w-b)(s-b) _ l 



wG 
1 



w 

4 x/w 



5(1 



1 



+ 2ssjxl-^ 
1 dT g _ C F a s (fi R )s 2 



(1 - G 

2 l-x 6 



1 + 6 

2w 



2^ 



G 1 — 2sx{, + 6 ^ G 



1 + 6 



r dx q 



7iwQG 
l-x„ 



(2s + 3w) 



4- 



X r. 



2w 
b 1 — x 5 



(l-x 6 )-l 

1 — X 

arcosh 



- 2sx g 



+ 



/3 2 (1 - 2sx 6 



(i-2sx g y 



(l-* 9 ) 2 



[s 2 (6s + w)x 2 g — 2sx g — wGq 



x 7sx„ — 8 H 



where 



r n 



G F m t 3 |V^ 



tb\ 



2ny/2 



-wQGq 



(9) 
(10) 



is the total decay width at LO and, in the notation of Ref. [9], 



m 



2 ■ 



s = hi+b-w), p = y^, 

2 s 



Q = 8y/l- ^ 



Gn 



1 + b - 2w + 



(1-6) 1 



10 



$(x 6 ) 



— (3 2 — artanh 



x b 



(11) 



Integrating dT b /dx b of Eq. (Q over x^, we recover the familiar result given by Eq. (2.8) 
in connection with Eqs. (2.2)-(2.4) and (2.6) of Ref. @j. 

As explained above, the GM-VFN results are obtained by matching the FFN ones to 
the ZM-VFN ones by subtraction, as 



1 dT,. 



r dx a 



GM-VFN 



Tq dx a 



1 dT r 



FFN 



where the subtraction terms are constructed as 



1 dT n 



1 (Tr. 



T Q dx a 



sub 



lim 

m6->0 \ Tq dx Q 



T dx a 



1 dT n 



sub 



FFN 



r dx a 



(12) 



(13) 



ZM-VFN 



W + -rest frame 



t-rest frame 



Figure 1: Definition of the polar angle 9 in the P^-boson rest frame. 



Taking the limit m b — > in Eq. (jUj), we recover Eq. (J7J) up to the terms 



Tq dxb 
1 rff, 



Yq d,Xn 



sub 



sub 



2tt 



C7 



1 + ^ 2 
1 - x b 



ln^|-21n(l-x 6 ) - 1 



m 



OL S (Hr) 1 + (1 - Xg) 



2tt 



In ^|-2 lux, - I ] . 

ml 



(14) 
(15) 



As already observed in Ref. [9], Eq. f fl4|) coincides with the perturbative FF of the tran- 
sition b -> b 1121. 



4 Analytic results for dT(t -t Bt+v t + X) 
I (dxB d cos 6) 

We now allow for the W boson in process to decay. For defmiteness, we consider its 
leptonic decay, which is cleaner than the hadronic one, so that we are dealing with the 
process 

t(Pt) -> b( Pb ) + W + (p w )(+g(p 9 )) B(p B ) +£+ + l y e + X, (16) 

for which we wish to calculate the doubly differential partial decay width, dY / (dxs d cos 9), 
at NLO. As mentioned above, 9 is the decay angle of the charged lepton £ + in the PF-boson 
rest frame (see Fig. CQ. 

As will be demonstrated in Sec. [5j the finite-m;, corrections are rather small in the 
case of process (JT]), much smaller than the contribution from gluon fragmentation. In 
the following, we, therefore, set m b = 0, i.e. we work in the ZM-VFN scheme, but we 
retain mg ^ and gluon fragmentation. Furthermore, we put p\ = m"^, i-e. we work 
in the narrow-width approximation. This allows us to employ the helicity density matrix 
formalism, so that the squared decay amplitude of process ffTB]) factorizes into the squared 
amplitude of process (pQ), with definite W-boson polarization, and the squared amplitude 
of the subsequent decay W + — > £ + + ve of the W boson just so polarized. The three 



8 



degrees of massive- vector-boson polarization, A = 0, ±1, are conveniently implemented by 
applying the covariant projection operators 



e^(±)e v *(±) 



w 




Pt - Pw 






Pw 




m 2 w 


1 

2 _ 


-sT 


1 PwPw 
m w 



Pw 



Pt 



PfPW y 
~2 PW 



m 



w 



w 

Pw 



ix _ PfPW fj, 
Pt _2 PW 



m 



w 



V _ PfPW y 

Pt ~Zj2 Pw 



IV 



m t \pw\ 



{Pt) P (pw)a 



(17) 



where e 0123 = 1 and \pw\ = \/ (Pt ■ Pw/ m t) 2 — my/ = m t\/W ~ s ( x b + x g )} 2 — w is the 
modulus of the VT-boson three-momentum in the top-quark rest frame. Performing the 
polarization sum, we recover the familiar completeness relation 



+i 



£y(A)e"*(A) = -cr + 



PwPw 



A=-l 



m 



w 



used in Sec. El 

Repeating the calculation of Sec. [3] at NLO in the ZM-VFN scheme with Eq. (JTTJ) 
instead of Eq. ffl8|) . we obtain 



d 2 T n 



dx a d cos 9 8 



-(1 + cosi 



)2^ + 3 (1 _ cos ^ )2 ^ 
dx n 8 dx„ 



- sm z 9 — -. 

4 dx n 



(19) 



where dT^/dx a are ^-independent coefficient functions listed in Eq. ( 1241) of the Appendix. 
As in Sec. [31 the top quark is assumed to be unpolarized. At LO, A = refers to 
the case when the top-quark spin is passed on to the bottom quark as is, while it is 
flipped for A = —1; A = +1 is prohibited by angular- momentum conservation, so that 
dT^/dXa vanishes at LO. At NLO, all values of A are allowed because of the presence of 
the additional spin-one boson g. 

There are two powerful checks for the correctness of Eq. (I24p . Firstly, integrating 
Eq. ffT9l) over cos#, we obtain 

(20) 



dY n 



dx a 



A=-l 



which agrees with Eq. ([7]) upon insertion of our expressions for dt^/dx a . Secondly, inte- 
grating dT^/dxb over Xb, we recover the results presented in Eqs. (15)— (17) of Ref. [2T] . 

The structure of the angular dependence of Eq. ffT9l) is preserved by the convolution 
with the FFs according to Eq. ([3D [or Eq. (|3J)], and we may project out the hadronic 
counterparts of dT^/dx a , 



dT x 

dxB 



E 

a=b,g Jl 



dx a dT X a (X B 
X a dX a V $a 



(21) 
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from the measured 9 dependence of d 2 



T/ (dxs dcosO) as 



i 



d 2 T 



d%B 



dcos9 f x (cos9) 



dxsdcosO 



(22) 



-i 



where 




- ±x + - 
2 2 



f°{x) = 2- 5x 2 . 



(23) 



In this way, we obtain three independent xb distributions, which we can use to constrain 
the 5-hadron FFs. 

5 Numerical analysis 

We are now in a position to explore the phenomenological consequences of our results 
by performing a numerical analysis. We adopt from Ref. [22J the input parameter values 
G F = 1.16637 x 1(T 5 GeV" 2 , m w = 80.399 GeV, m t = 172.9 GeV, m b = 4.78 GeV, 
m B = 5.279 GeV, and \V tb \ = 0.999152. We evaluate u^Qir) at NLO in the MS scheme 
using Eq. (9.5) of Ref. [22] , retaining only the first two terms within the parentheses, with 

(5) 

rif — 5 active quark flavors and asymptotic scale parameter A^ = 231 MeV adjusted 

such that ots\mz) = 0.1184 for rriz = 91.1876 GeV [22]. We employ the nonperturbative 
S-hadron FFs that were determined at NLO in the ZM-VFN scheme through a joint fit 
[IS] to e + e"-annihilation data taken by ALEPH [23] and OPAL [24] at CERN LEP1 and 
by SLD [25J at SLAC SLC. Specifically, the power ansatz D^z,^ 1 ™) = Nz a (l — z) 13 was 
used as the initial condition for the b — >■ B FF at /i™ 1 = 4.5 GeV, while the gluon and 
light-quark FFs were generated via the DGLAP evolution. The fit yielded N = 4684.1, 
a = 16.87, and (3 = 2.628 with Xd.o.f. = 1-495. We choose fiR = [If = rn t . 

We first consider the quantity dT(t — > B + X)/dxs taking the W boson to be stable. 
Our most reliable prediction for it is made at NLO in the GM-VFN scheme and includes 
finite-m^ corrections. For the time being, we implement the latter in terms of the energy 
scaling variable z = E B /E a using Eq. ([3]). In Fig. [2J we study for this prediction the size 
of the NLO corrections, by comparing the LO (dotted line) and NLO (solid line) results, 
and the relative importance of the b — >■ B (dashed line) and g — )■ B (dot-dashed line) 
fragmentation channels at NLO. In order to expose the size of the NLO corrections at 
the parton level, we evaluate the LO result using the same NLO FFs. We observe from 
Fig. [2] that the NLO corrections lead to a significant enhancement of the partial decay 
width in the peak region and above, by as much as 25%, at the expense of a depletion 
in the lower-x^ range. Furthermore, the peak position is shifted towards higher values of 
xr- The g — > B contribution is throughout negative and appreciable only in the low-a;^ 
region, for xb^0.2. For higher values of xb the NLO result is practically exhausted by 
the b — > B contribution, as expected [9]. In the following, we only consider NLO results 
unless otherwise stated. 

In Fig. [31 we study the improvement over previous calculations, at NLO in the ZM- 
VFN scheme neglecting g — >■ B fragmentation and finite-m^ corrections (dotted line), 
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Figure 2: dT(t — > B + X)/dxB as a function of xb in the GM-VFN (m& 7^ 0) scheme 
including finite-m^ corrections implemented with the energy scaling variable z = Es/Ea 
using Eq. Q. The NLO result (solid line) is compared to the LO one (dotted line) and 
broken up into the contributions due to b — > B (dashed line) and g —> B (dot-dashed 
line) fragmentation. 
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X B 

Figure 3: dT(t B + X)/dx B as a function of x B at NLO. The GM-VFN (m b ^ 0) result 
including finite-m^ corrections implemented with the energy scaling variable z = E B /E a 
using Eq. ([3]) (solid line) is compared to the ZM-VFN (m& = 0) results for m B = 
excluding (dotted line) and including (dashed line) g —> B fragmentation. 



gained by switching to our preferred mode of evaluation, at NLO in the GM-VFN scheme 
including g — > B fragmentation and finite-m^ corrections (solid line). For comparison, 
we also show the full ZM-VFN prediction, including g — > B fragmentation, for = 
(dashed line). We observe from Fig. [3]that the improvement is threefold. First, the finite- 
m B corrections are responsible for the appearance of the threshold at x B = Pb = 0.07783. 
Second, g — > B fragmentation leads to a significant reduction in size in the threshold 
region, as is familiar from Fig. [2j Third, the finite-m& corrections lead to a moderate 
reduction in size throughout the whole x B range allowed. A similar observation was made 
for the Xb distribution in the perturbative-FF approach [9]. 

For a more quantitative interpretation of Fig. [31 we consider in Fig. H] the ZM-VFN 
result for m B = without g — )■ B fragmentation (dotted line) and the full GM-VFN 
result implemented with the energy scaling variable (solid line), both normalized to the 
full ZM-VFN result for m B = 0. We observe that the omission of g — > B fragmentation 
entails an excess by a factor of up to 2 close to threshold, while the finite-m 6 corrections 
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Figure 4: dT(t B + X)/dx B as a function of x B at NLO in the GM-VFN (m b ^ 0) 
scheme including finite-m^ corrections. The implementation with the energy scaling 
variable z = E B /E a via Eq. (J3J) (solid line) is compared to the one with the light-cone- 
momentum scaling variable z — VbIPo, v ^ a Eq- © (dashed line). For comparison, also the 
ZM-VFN (rrib = 0) result for mg = excluding g — > B fragmentation is shown (dotted 
line). All results are normalized to the ZM-VFN (m^ = 0) one for tub = including 
g — > B fragmentation. 



just amount to a few percent. This motivates us to adopt the ZM-VFN scheme including 
g —7- B fragmentation in the treatment of d 2 T{t — > B£ + i>z + X) / (d%B dcosO) below. 

Figure 0] also includes the result obtained in the GM-VFN scheme implemented with 
the light-cone-momentum scaling variable z = VbIvX via Eq. fllj (dashed line). Compari- 
son with its counterpart for the energy scaling variable allows us to assess the theoretical 
uncertainty in the implementation of finite-m^ corrections in the GM-VFN scheme. We 
observe from Fig. H] that the two implementations of finite-m^ corrections lead to very 
similar results over most of the xb range. Appreciable differences only appear close to the 
kinematical bounds. However, these should be taken with a grain of salt because the FFs 
used here [18] refer to the ZM-VFN scheme and do not carry any specific information on 
finite-ms effects and the scheme of their implementation. 

In the remainder of this section, we consider the quantity d 2 T(t — >■ B£ + i>e + X) 
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X B 



Figure 5: dT x (t — > Bl^vi + X)/dxB as functions of xb in the ZM-VFN (nib = 0) scheme 
including finite-m^ corrections implemented with the energy scaling variable z = EB/E a 
using Eq. ([3]). The NLO results for A = (solid line), A = — 1 (dashed line), and A = +1 
(dot-dashed line) and the LO results for A = (upper dotted line) and A = — 1 (lower 
dotted line) are compared to each other. 



jidxB dcosd), which we study at NLO in the ZM-VFN scheme including g — > B fragmen- 
tation and finite-m^ corrections evaluated using the energy scaling variable. In Fig. |5l 
we present the xb distributions of dr x /dxs for A = (solid line), A = — 1 (dashed line), 
and A = +1 (dot-dashed line). For comparison, also the LO results for A = 0, —1 (dotted 
lines) are shown. As may be gleaned from Eq. (1241) . the latter amount to respectively 
1/(1 + 2w) = 69.59% and 2w/(l + 2w) = 30.41% of the LO result. In these two cases, 
the NLO corrections are similar in size and shape to the unpolarized case (see Fig. 
As explained in Sec. HI A = +1 is prohibited at LO, which explains the smallness of the 
corresponding result. 
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6 Conclusions 



Let alone its discovery potential with regard to the Higgs boson and new physics beyond 
the standard model, the LHC is a formidable top factory, which, among other things, 
will allow for the study of the dominant decay mode t — > BW + + X with unprecedented 
precision in the long run. In particular, this will enable us to deepen our understanding of 
the nonperturbative aspects of B-hadron formation by hadronization and to pin down the 
b — >■ B and g — >■ B FFs. The key observable for this purpose is the xb distribution dT/dxB 
of t — > B + X. By measuring the angular distribution of the VT-boson decay products, 
also the three components dT x /dxB of dT/dxs corresponding to the polarization states 
A = 0, ±1 of the W boson may be determined, which constrain these FFs even further. 

We studied the quantity dT/dx B at NLO in the GM-VFN scheme [13115116117115] 
using nonperturbative 5-hadron FFs determined by a global fit [18] of experimental data 
from the Z factories [23. 2 ^125] . relying on their universality and scaling violations [T9] . 
This allowed us to investigate for the first time finite-mb corrections and the contribution 
from gluon fragmentation to dr/dxs- We also analyzed the size of finite-m^ effects and 
their theoretical uncertainty due to the freedom in the choice of the scaling variable. Since 
the finite-m;, effects turned out to so moderate, we neglected them, for simplicity, in our 
study of dr x /dxs, which we treated at NLO in the ZM-VFN scheme taking into account 
finite-m^ effects. 

Comparing future measurements of dT/dxs and dT x /dxB at the LHC with our NLO 
predictions, we will be able to test the universality and scaling violations of the 5-hadron 
FFs. These measurements will ultimately be the primary source of information on the 
.B-hadron FFs. The formalism elaborated here is also applicable to the production of 
hadron species other than B hadrons, e.g. pions, kaons, protons, D mesons, etc., through 
top-quark decay. 
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Appendix 

The coefficient functions f£ in Eq. ( Fl9l) exhibit the structure 
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where 



Pqq{ x ) — Cp 



1 + x 2 
1 — x 
l + (l-x) 2 
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(25) 



are the timelike q — > q and q g splitting functions at LO, respectively. Introducing the 
short-hand notation 

w± = \± y/w, (26) 



and defining 
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(27) 
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where 9(x) = J^dtS^t) is the Heaviside step function, we have 
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